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NOMENCLATURE 


Cn(S)  n  times  continuously  differentiable  on  S 

p  pressure 

U  component  velocity  in  x-direction 

v  component  velocity  in  y-direction 

M  1  32 [  ]/3x2  +  32 [  ]/32[  ]/3y2 

AM  ]  3“[  ]/3x2  +  2  5!|iLl  +  ]/3y“ 

3x  3y 

3S  boundary  of  domain  S 

n.  £•  u 

■»  implies 
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/  SUMMARY 

' 

r  For  more  than  the  last  20  years  there  has  been  a  concerted  effort  to 
solve  the  stationary  Navier-Stokes  equations;  however,  this  has  only  been  suc¬ 
cessful  for  a  few  special  cases  of  primarily  academic  interest.  An  alternative 
approach  has  been  to  solve  the  equations  numerically,  and  then  compare  the 
results  with  experiment.  On  occasion,  such  comparisons  are  in  good  agreement. 
However,  such  results  are  of  dubious  value  since  one  has  no  a-priori  way  of 
knowing  the  relevance  of  such  results  until  they  are  explicitly  compared 
against  experiment.  Therefore,  it  would  seem  reasonable  to  conclude  that  the 
present  approaches  to  solving  the  Navier-Stokes  equations  are  of  limited 
value.  Accordingly,  it  is  the  purpose  of  this  paper  to  show  that  there  does, 
indeed,  exist  an  equivalent  representation  of  the  problem  that  has  significant 
potential  in  solving  such  problems.  This  is  due  to  the  fact  that  this  equiva¬ 
lent  representation  of  the  problem  consists  of  a  sequence  of  Fredholm  Integral 
Equations  of  the  second  kind,  and  the  solving  of  this  type  of  problem  is  very 
well  developed.  In  addition,  for  the  problem  in  this  form,  there  is  an  excel¬ 
lent  chance  to  also  determine  explicit  error  estimates,  since  one  would  now  be 
dealing  with  bounded  linear  operators,  rather  than  unbounded.  s 


INTRODUCTION 


Ideally,  one  would  like  to  obtain  a  numerical  solution  to  the  Stationary 
Navier-Stokes  problem  that  is  within  some  prescribed  degree  of  accuracy  of  the 
true  solution.  Unfortunately,  it  is  not  possible  to  accomplish  this  objective 
by  just  applying  some  existing  numerical  technique  directly  to  the  Stationary 
Navier-Stokes  equations,  since  error  analysis  for  nonlinear  equations  is,  for 
all  practical  purposes,  nonexistent.  However,  since  it  is  possible  to  replace 
the  Stationary  Navier-Stokes  problem  by  an  equivalent  sequence  of  linear 
partial  differential  equations  (ref.  1),  providing  some  rather  general  condi¬ 
tions  can  be  met,  it  is  at  least  theoretically  possible  to  obtain  error  esti¬ 
mates  due  to  the  vast  wealth  of  knowledge  known  about  linear  equations.  How¬ 
ever,  even  though  attaining  such  error  estimates  is  within  the  realm  of 
possibility,  one  should  not  be  misled  into  thinking  that  such  a  task  is  a 
small  undertaking,  for  this  is  most  certainly  not  the  case.  However,  the 


i 


chances  of  someone  accomplishing  this  task  are  significantly  greater  if  the 
basic  sequential  problem  is  transformed  into  an  equivalent  form  that  is  more 
suitable  to  error  analysis.  Accordingly,  in  this  paper  it  will  be  proven  that 
the  said  sequential  problem  can  be  transformed  into  a  form  that  has  signifi¬ 
cantly  greater  changes  of  yielding  explicit  error  estimates. 

GENERAL  DEVELOPMENT 


(1) 


(2) 


Mas)  *  t>l  ,  M3s>  -  t>2 

*  y 

Let 

vAAO  +  0  AO  -  0  +  f,  ~  f2 

y  a  x  y  y  x 

be  denoted  by  P(0)-  Hence  (2)  can  be  expressed  by 

P(o)  -  o  ,  0X(3S)  «  b,  ,  ♦  OS)  *  b2 

where  P  can  be  interpreted  as  a  mapping  from  CI,(S)  into  C°(S).  As  demon¬ 
strated  in  reference  1,  solving  (2)  is  equivalent  to  solving  the  sequence  of 
equations 

P(0  )  +  P* (0 
n  n 

3_K 

3X  dS  ' 


)  (il'  —  )  *  0 


»,  •  vl„  •  b‘  for  "  ■  •  • 


The  Navier-Stokes  equations  are  given  by 

uu_.  +  vu  +  p  -  v  Au  +  f1(x,y)  «  0  \ 

y  |  in  S 

uvx  +  wy  +  py  -  v  Av  +  f 2  (x,y)  *  0  J 

with  boundary  conditions 

u(3S)  «  -b 2  (3S)  ,  v(3S)  =  bx  , 

where  S  is  a  two  dimensional  Green's  domain  with  surface  3S, 

f i (x,y)  e  C1(S)  ,  and  f2(x,y)  e  C^S)  , 

or  equivalently, 

vAAO  +  0V  ~  0V  AO  +  f.  -  f,  *  0  in  S 

y  a  Ay  x 
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providing  certain  conditions  are  satisfied  (ref.  1).  (i.e.,  primarily  that 
t|»0,  the  initial  guess  is  reasonably  good),  for  further  Insight  into  the  sig¬ 
nificant  latitude  on  ip0,  the  reader  should  consult  reference  2. 


P(tf»0)i  ]  -  vAA(  ]+!/<„  A[  ]x  +  A4>6  f  ]  -  A*0  (  ]x  -  A(  ] 


Equivalent  to  equation  (3)  is  the  problem 


P<V  +  p,<*n>*n-°  *  TT 


as 


o  , 


3^ 

3y 


•  0  ,  n  *  0,1,2, 
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by  letting 
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n+i 


Now  consider  (4)  under  a  slight  change  in  boundary  conditions,  in  particular: 


*(♦_ )  +  P *  (*„)**  '  0 

m  mm 


_ ffl 

3n 


as 


oi  *  0|1|2)  •  •  • 


(5) 


If  is  a  solution  to  (5),  it  directly  follows  that  is  also  a  solution 

of  (4),  and  vice  versa  (i.e.,  equation  (4)  ■*  ym|  *  constant,  without  loss  of 

I  3S 

generality  assume  constant  *  0).  Therefore,  it  is  sufficient  to  focus  our 
attention  on  equation  (5),  which  in  detail  is  given  by 


uAAiJ/  + 
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m.. 


SS 


-  *m  Mm  +  P(*m>  •  0 
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With 


*nOS)  *  0 


0  ,  m  -  0,1,2,  .  .  . 


Hence, 

"+m  *  v  J[  G(\  A\  +  '  \A\)ds  +ifsPG  ds  (6) 

where  G  is  the  Green’s  function  of  the  biharmonic  Equation  [3] 

AA$  *  0  in  S,  with  $(3S)  *  0  ,  $  (3S)  ■  0 

n 

Lemma: 


Proof. O 


(6)  -  (5) 


For  convenience  and  clarity,  equation  (6)  will  be  expressed  in  the  form: 

$  •  A(4in)^  +  f<*m)  (6) 

where 

*<*„>'  1  -  i  f  0(*  A|  ]x  t  4*^1  )y  -  A*.  I  )„  -  *.  A[  )y)d» 

and 

f(*m)  -  -  J  PU>JG  ds 
m  v  's  m 

Lemma:  If  the  hypothesis  of  theorem  1  of  reference  1  is  satisfied,  then  there 
exists  a  solution  of  equation  (6)  for  all  m  and 


n 

4>*  -  Urn  £  -  *# 

n-»®  i*o 

where  ij<*  is  the  unique  of  equation  (2). 

Theorem:  Under  the  condition  of  the  above  lemma,  is  a  solution  to  the 
Fredholm  Integral  Equation  of  the  second  kind 


^*(x',y')  +  J  ^(x.yJK^x*  ,y‘  ,x,y)dx  dy  «=  fU^Hx'.y') 
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where 
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-AGy (x' ,y ' 


»x>yHm  <x*y>  +  4GxV 
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Proof:  Since  will  be  a  solution  of  (6)  under  the  hypothesis,  it  is  suffi¬ 

cient  to  show  that  a  solution  to  (6)  is  alos  a  solution  of  (7).  For  clarity, 
denote  (6)  in  the  form 


-  A(^0H*  *=  f(<{>0) 


or  equivalently. 
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*  _ 


o  **  " 
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By  Green’s  Lemma 
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-a(vO  ■  -(v.,)xx  -  (s*.x)„ 

-  ■  cy*!xx\  +  ('Cyy!=x  ■  cy*‘x,\ 
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CONCLUSION 

Therefore,  under  the  conditions  as  cited,  the  sequence  of  solutions 
generated  by  (7)  converge  to  .the  solution  of  (1)  in  the  following  sense 

"t  », 

i°o 

fl  *  if>*  -  n, 

o  A 

Therefore, 


■  *0  +  £ 


Hence  the  original  problem,  equation  (1),  which  falls  within  the  frame¬ 
work  of  nonlinear  operator  theory,  an  area  that  little  is  known,  has  been 
replaced  by  equation  (5),  which  falls  within  the  framework  of  bounded  linear 
operator  theory,  an  area  for  which  there  exists  a  vast  wealth  of  information. 

Of  course,  it  could  be  argued  that  it  is  theoretically  possible  that  the 
sequential  representation  of  the  solution  converges  so  slowly  that  the  results 
are  of  questionable  value.  However,  from  all  indications  this  will  not  be  a 
problem,  for  as  pointed  out  in  reference  1,  there  are  parameters  at  one's 
disposal  in  the  method  that  can  be  adjusted  to  speed  up  the  rate  of  conver¬ 
gence;  in  fact  it  was  demonstrated  in  reference  2  that  by  judicious  selection 
of  the  variable  parameters  available  they  were  able  to  get  numerically  adequate 
convergence  with  just  a  few  iterations. 
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